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, In this talk I describe some applications of random matrix models 

0^ ' to the study of A/" = 1 supersymmetric Yang-Mills theories with matter 

I fields in the fundamental representation. I review the derivation of the 

• Veneziano-Yankielowicz-Taylor/Affleck-Dine-Seiberg superpotentials from 
constrained random matrix models (hep-th/0211082), a field theoretical 
justification of the logarithmic matter contribution to the Veneziano-Yankie- 
lowicz-Taylor superpotential (hep-th/0306242) and the random matrix based 

Oh, solution of the complete factorization problem of Scibcrg-Witten curves for 

^ • N = 2 theories with fundamental matter (hep-th/0212212). 



1. Introduction 

One of the most important and challenging problems in physics is the 
understanding of the non-perturbative properties of gauge theories. In the 
last 10 years a series of significant breakthroughs were made in the study 
of supersymmetric gauge theories, where the additional symmetry proper- 
ties allowed to obtain exact non-perturbative information at the same time 
unraveling various unexpected mathematical structures. 

The first wave of research was based on exploiting the holomorphic prop- 
erties of various J\f = 1 non-perturbative quantities (see e.g. [1]). The sec- 
ond breakthrough, this time giving very complete information on A/" = 2 
theories, showed that the low energy dynamics is encoded in properties of 
certain elliptic (or hyperelliptic) curves — the Seiberg-Witten curves [2, 3]. 

The most recent breakthrough arose in the wake of the gauge the- 
ory/string correspondence [4]. Using D-brane constructions of gauge theo- 
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ries and performing 'geometric transitions' to a dual geometric setup, super- 
potentials in supcrsymmetric gauge theories where related to the geometry 
of the dual Calabi-YaTi manifold [5, 6]. In a subsequent crucial development 
Dijkgraaf and Vafa interpreted the geometrical formulas in terms of a saddle 
point solution of an associated random matrix model [7, 8]. Subsequently 
the resulting random matrix prescription for calculating the superpotential 
was proven purely in the field-theoretical context by using diagrammatic 
analysis [9] and anomaly considerations [10]. 

In this talk, after recalling some basic facts about Af = 1 SYM theories, 
I would like to review three developments [11, 12, 13] in the extension of 
the original Dijkgraaf- Vafa proposal to supcrsymmetric gauge theories with 
matter fields in the fundamental representation. 

The original proposal involved only superpotentials expressed in terms 
of the glucball (gaugino bilinear) superfield S. In [11] wc showed how one 
could incorporate into the framework mesonic superfields, and how to obtain 
the classical Veneziano-Yankielowicz-Taylor/Affleck-Dine-Seiberg superpo- 
tentials directly from the matrix model. 

The structure of the superpotentials involve typically logarithmic and 
power-series terms. The latter are quite well understood [9, 10] while the 
former remain somewhat mysterious. In [12] we gave a field-theoretic di- 
agrammatic derivation of the matter induced logarithmic part of the VYT 
potential. 

Finally I would like to review the use of random matrix models to obtain 

an explicit complete factorization [13] of the Seiberg-Wittcn curve of an 
J\f = 2 U{Nc) SYM theory with Nf < Nc fundamental flavours. This 
is equivalent to finding the submanifold of A/" = 2 moduli space of vacua 
where all monopoles in the theory become massless. 



2. Physics of A/" = 1 SYM theories with fundamental matter 

In order to study the dynamics of the theory at low energies one is in- 
terested in obtaining the low energy effective action for the relevant degrees 
of freedom. In the case of = 1 gauge theories most of such degrees of 
freedom may be described by chiral (and anti-chiral) superfields. Then the 
constraints of super symmetry restrict the effective action to be of the form: 

Seff = j d^6d^eK{^,^)+ J d'^eWeff{^) + J d^eWeffi^) (1) 

where $ stands here for a generic chiral superfield. In general little is 
known about the non-chiral part but a lot of information can be 

obtained on the effective superpotential We//(^)- Its knowledge allows us 
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to determine vacuum expectation values (VEV) of $ from the equation 

dWeffi^) 



(2) 



Let us briefly recall what is known about effective supcrpotcntials first for 
pure M = 1 SYM with gauge group SU{Nc), and then for the theory with 
Nf fundamental flavours. 

Pure ^^ =1 SYM 

A natural chiral superfield which can be built from the gauge field is 

>5=-^trWaW" (3) 

whose lowest component is the gaugino bilinear. The form of the effective 
superpotential Weff{S) has been determined by Veneziano and Yankielow- 
icz based on anomaly considerations [14]: 

Wvy{S) = -S log ^ + 2TriToS (4) 

The second piece is just the tree level YM coupling (tq = ^/27r + ATri/gyM) 
and can be absorbed into the first one through a redefinition of A. The 
power S^" is determined by the anomaly, while the exponent of A is just the 
coefficient of the 1-loop /3 function, in order for the superpotential to be RG 
invariant. Wvy{S) leads to a nonzero VEV for S — gaugino condensation. 

Af — 1 SYM with Nf < JVc flavours 

When we add to the theory matter fields in the fundamental represen- 
tation there appear additional mesonic chiral superfields: 

^ij = QiQj (5) 

Veneziano, Yankielowicz and Taylor determined the appropriate superpo- 
tential again using anomaly considerations [15] 

^^ /.r S ^, detAT ^, fS^--^f detX\ 

Wvyt{S,X) = {Nf-N,)S\ogj^-S\og^^ = -Slog ^3^^_^^ J 

(6) 

Note the modification of the anomaly coefficient {Nc — Nf) and the 1-loop 
P function 3Nc — Nf due to the matter fields. Later we will show how these 
modifications arise from the matrix model framework. After one integrates 
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out S one is left with the Affleck-Dine-Seiberg superpotential [16] for the 
mesonic superfield only: 

W,ns(X) = (N,-N,){^-^) (7) 

Intrilligator-Leigh-Seiberg (ILS) linecirity principle 

Suppose that the gauge theory is deformed by some tree level poten- 
tial Vtree{X) (c.g. mass terms for the flavours and some self-interactions). 
Intrilligator, Leigh and Seiberg [17] conjectured that the full dynamics is 
described just by the sum 

WvYT{S,X)+Vtree{X) (8) 

i.e. the dynamical potential is not modified by the deformation. 

This concludes the very brief review of classical results on jV = 1 SYM 
theories. In the next section we will present the main new ingredient — 
the link with random matrix models proposed by Dijkgraaf and Vafa, and 
proceed to describe its applications. 



3. The Dijkgraaf- Vafa correspondence 

Dijkgraaf and Vafa proposed an effective way of calculating superpo- 
tentials for A/" = 1 theories with adjoint fields and arbitrary tree level su- 
perpotentials [8] (later generalized also to include fields in the fundamental 
representation [18]): 

W,ff{S) = Wvy{S) + iVc^^^^ + ^x=i(^) (9) 
where the J^iS are defined through a matrix integral 



e 



J D^DQiDQiexp!^-jWtree{^,Qi,Qi)'^ (10) 



Here $ is an TV x TV matrix, while the QiS are N component (complex) 
vectors. In this expression one takes the limit oo, gs ^ with 

QsN = S = const in order to isolate graphs with the topology of a sphere 
{Ty^=2) and of a disk {T^=i). 

The above expression gives a prescription for the superpotential only in 
terms of the glueball superfield S. However it is interesting to ask if one 
could obtain directly from matrix models the superpotential which involves 
also mesonic superfields. A proposal for doing that was given in [11] and 
will be described in the next section. 
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4. Mesonic superpotentials from Wishart random matrices 

In order to express the effective superpotential in terms of mesonic su- 
perfields Xij = QIQj, it was proposed in [11] to perform only a partial 

integration over the Q's in (10) and impose the constraint Xij = qIQj 
directly in the matrix integral i.e. 

e-^AW^x{S,x) ^ J DQiDQi6iXij-QlQj)exp!^-^Vtree{X)'^ (11) 

Then the effective superpotential involving both S and X is obtained from 

WeffiS, X) = Wvy{S) + Nc ^^^'^lf^^ + ^x=i{S, X) (12) 

For theories with matter fields only in the fundamental representation (as 
written in (11)) the above simplifies since Tx='^ ~ ^■ 

Prom (11) we see that Vtreei^) contributes directly to J^y^=i{S,X). This 
is in complete agreement with the ILS principle where the 'dynamical' 
Veneziano-Yankielowicz-Taylor superpotential is not influenced by the tree 
level deformation. 

The 'dynamical' contribution to J^x=i{S,X) will come from the con- 
strained integral over Nf vectors of length N 



I 



DQiDQiS{Xij-QlQj) (13) 



Up to an inessential term exp(— tr X) this is just the probablity distribution 
of (complex) Wishart random matrices [19]. The result is known^ and for 
the case N > Nf reads 



/ 



(2vr) 



N(N+1) 
2 



DQDQSiQQ -X)= ' > (detX)^^"^^/ (14) 

ni=Ar-JV^+lU -1)' 

Extracting now the large A'^ asymptotics according to (11) one sees that 
the normalization factor gives rise to a NfS log S term, which reflects the 
contribution of the matter fields to the anomaly. The detX term leads to 
the correct dependence on the mesonic superfield. The full superpotential 
is then 

WvYiS)+NfSlog ^-51og +Vtree{X) = WvYT + VtreeiX) (15) 

as expected from anomaly considerations [15]. It is quite surprising that 
both terms arise from the classical random matrix Wishart distribution. 



^ See e.g. [20] for a general proof and eq (15) in [21] for the numerical coefficient. 
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The case with Nf = Nc 

It is interesting to consider the case Nf = Nc with Vtree{^) = 0. From 
the general form of (15) we see that the logarithmic term vanishes and S 
appears only linearly as 

Weff = S log (det X/A'^^-) (16) 

thus it generates in a natural way a constraint surface (moduli space) sat- 
isfying the equation detX = A^^". On this surface the effective poten- 
tial vanishes. This is the correct behaviour for Nf = Nc and agrees with 
Seiberg's quantum constraint 

det X - BB = A'^^'^ (17) 

when the baryonic fields are integrated out [22]. However a direct inclu- 
sion of the baryonic superfields into the random matrix or a more general 
combinatorial framework remains an open problem [23, 24]. 



5. Perturbative justification of mesonic superpotentials 

For adjoint matter fields the non-logarithmic part of the Dijkgraaf-Vafa 
prescription has been verified by a perturbative diagrammatic calculation 
[9] and arguments based on generalized Konishi anomalies [10]. However 
the origin of the logarithmic Veneziano-Yankielowicz term remains from 
this point of view quite obscure — as it should arise from Feynman graphs 
which involve 'dynamical' vector superfield loops. 

Since a logarithmic term also appeared from a random matrix calcula- 
tion in the constrained matrix integral it is interesting to check if the matter 
induced part of the Veneziano-Yankielowicz superpotential may be also ob- 
tained diagramatically directly in the gauge theory. This calculation was 
done in [12]. 

The matter contribution to the effective potential in S (with all matter 
superfields integrated out) is given by [9] 

J DQDQ '^"'"^^'^ {-^Q{a-iW°'dc)Q+Wtree{Q,Q)) j^^gj 

where W" is the external field related to S through (3) and da = gf^- In 
[12], in order to obtain the effective superpotential involving also mesonic 
superfields we introduce the superspace constraint 

X = QQ (19) 

by inserting into (18) a Lagrange multiplier chiral superfield a. Using the 
fact, exploited in [9], that the antichiral sector does not influence chiral 
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Fig. 1. Only tree level graphs survive (like the one shown in fig. Ic). 
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Fig. 2. The Schwinger-Dyson equation for F. 

superpotentials we introduce also an antichiral partner a with a tree level 
potential i.e. 

J ct^xd^e aa + J ct^xSO c? . (20) 

The path integral over a. can be carried out exactly and yields (c.f. [9]) 

J d^x(fe aOa. (21) 

Therefore the final path integral which one has to evaluate is 

J DaDQDQ (22) 

This is a nontrivial interacting field theory, but since we want to extract 
only the trW^ terms we can allow only at most two W insertions in a 
QQ loop. The structure of the integration over fermionic momenta then 
significantly reduces the number of contributing graphs. In particular we 
are left only with graphs coming from (22) which have the structure of QQ 
loops connected by at most one a propagator and a propagators connected 
to the external field X as shown in fig. 1. Since the a propagators are then 
evaluated necessarily at zero momentum one has to include an IR cut-off 

The form of the superpotential then follows from (for details see [12]) 

Slogdet^ (23) 

where F is the a 1-point function determined by a Schwinger-Dyson equa- 
tion (fig. 2): 
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The only IR-finite solution is F = SX ^ which when inserted into (23) 
yields the expected result (compare with (15)) 

NfS log ^ - S log det ^ (25) 
including the matter generated logarithmic term. 



6. Random matrices and Seiberg-Witten curves 

In the preceeding sections we described the rcdcrivation using either 
random matrix or diagrammatic methods of the Veneziano-Yankielowicz- 
Taylor/AfHeck-Dine-Seiberg superpotentials. It is interesting to ask whether 
one can use the close relation between random matrix models and SYM 
theories with matter fields to obtain some new nontrivial information about 
the latter. An example of this kind of result is the effective factorization of 
Seiberg-Witten curves from random matrix considerations [13]. Before we 
describe it in the following section let us first recall some basic facts about 
J\f = 2 SYM theories. 

At low energies M =2 U{Nc) SYM theories develop an iVc-dimensional 
moduli space of vacua parametrized by vacuum expectation values 

«P=(^tr$^^ (26) 

with 1 < p < A'^c- The IR dynamics of the theory around each such vacuum 
is described by the Seiberg-Witten curve 

y2 = P^^{x, ukf - 4A2^-^/ Hix + rrii) (27) 

i=l 

where the polynomial PNc{x,Uk) = {det{xl — ^)) = X]a=o '^a^^''"" 
pends in an explicit way on the Uk's: 

a 

asa + ^ kSa-kUk = (28) 

k=0 

so = 1, uo = (29) 

Of particular interest (cf. [25]) are the vacua where all the monopoles of the 
theory are massless. Then the Seiberg-Witten curve completely factorizes, 
i.e. the right hand side of (27) has only two single zeroes, all the remaining 
zeroes are double. The Seiberg-Witten curve then has to have the form 



= (x - a){x - b)HN^^i{xf 



(30) 
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The problem is to find the moduh {ti^j's for which this happens. One 
expects a 1-parameter family of solutions. For pure Af = 2 theory {Nj = 0) 
the explicit form of this submanifold of vacua has been found using special 
properties of Chebyshev polynomials [25]. Unfortunately these methods 
cannot be extended to the case of Nf > 0. In [13], using random matrix 
model techniques, we obtained formulas for when the general form of the 
Seiberg-Witten curve (27) completely factorizes for Nf < Nc and arbitrary 
masses m^. 

We note that the general problem of factorizing (27) is highly nonlinear 
and involves coupled sets of polynomial equations. In fact we would not 
expect a-priori that a closed-form analytic solution would exist at all. 

In the following section we will briefly describe the results of [13]. Other 
studies of the close link between Seiberg-Witen curves and random matrices 
include [26, 27, 28, 29, 30]. 



7. Factorization solution of Seiberg-Witten curves for theories 
with fundamental matter 

The main tool which allows to obtain the factorization of SW curves is 
the study of the M = 2 theory deformed by a tree-level superpotential 



We 

Wtree = Vffp--tr$f. (31) 
^ P 

Then once the factorization solution Up"^*' is known, the effective superpo- 
tential is given by 

Nc 

Weff = J2 m, T), (32) 

p=l 

which should then be minimized with respect to the parameter T of Up"''^' . 
If wc integrate in S by performing a Legendre transformation with respect 
to log A^^'^"'^-'" we obtain the superpotential 

H^e//(-5,txi,J^,A)=51ogA2^=-^/ + We//(5,txi,i^) = 

Nc 

=Slog A2^=-^/ - S\og^'^^^-^f + Y^gpuf/^*-{VL,ui){2,2.) 

p=i 

In order to obtain the factorization solution Up"^*' from random matrix 
models we should first use the Dijkgraaf-Vafa prescription to obtain Weff{S) 
from a random matrix expression and then recast it in the gauge-theoretic 
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form Weff{S,ui,il) defined by (33) which is linear in the couplings gp — 
this is the most difficult part of the computation. Then one can directly 
read off the factorization solution from the coefficients of the couplings. 
The structure of the solution obtained in [13] is the following: 



^fact. ^ NM^'^iR, T)+Y, Wp"^''""'(i?, T, rrii) (34) 



i=l 



where the two random matrix parameters R and T are related to the physical 
parameters A, ui by the constraints 

A^^-^/ = 7 ^ (35) 



N; 

m = A/"cr - ^^mi + T- V(wi + r)2 -4i? (36) 

In (34) Up^^^{R, T) is the factorization solution for pure M = 2 theory 

1 b/2] . \ \ 
U^^%R, = - E ( ( i?9rP-2? (37) 

while the explicit form oi U^"'^^'^^ {R, T, rrii) is given in section 7 of [13]. Here 
we just cite the result for p < 3: 

^rnatter = ^ _ + ^(^ + T)"^ " (38) 

Umatter^^^^ T, m) = ^ (n? - 2R - + (T - m)v'(m + T)2 -4i?j (39) 



W^""^^(i?, T,m) = l{ -m^ - 6RT - + (m^ + 2R-mT + T^) 
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(40) 



Let us note some striking features of the solution (34). Firstly, it has 
an extremely simple dependence on the number of colours. Secondly, each 
flavour contributes linearly. Thirdly, the solution for pure SYM theory 
appears as a part of the expression. All of these features are quite surprising 
if we keep in mind the very much nonlinear character of the factorization 
problem for the curve (27). Also from the physical point of view such 
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unexpected linearization seems to suggest some hidden structure of the SYM 
theory with fundamental flavours. 

Recently, the factorization solution (34) was used in [31] to rederive from 
the Seiberg-Witten curve perspective the Afileck-Dine-Seiberg superpoten- 
tial (7). 

8. Discussion 

The link of random matrices to supersymmetric gauge theories is one of 
the most exciting theoretical developments in the last year. It is notable as 
giving a completely new exact application of random matrix models. The in- 
terest lies both in reinterpreting classical random matrix results/ensembles 
in a new langage and context, and in Tising the novel random matrix meth- 
ods to obtain new results in supersymmetric gauge theories. It raises also 
numerous questions of a more mathematical nature. In particular it would 
be interesting to understand the precise interrelation between Calabi-Yau 
constructions and random matrices in the most general setting. This link 
was in fact at the origin of the Dijkgraaf-Vafa proposal. 

Last but not least it would be interesting to understand the mathemati- 
cal structures which allowed to use random matrix model calculations to fac- 
torize Seiberg-Witen curves. The solution obtained in [13] is based on a rein- 
terpretation of a random matrix formula using input from gauge-theoretical 
reasonings involving such non-perturbative concepts as 'integrating-in' etc. . 
It would be fascinating to link directly the underlying geometry of factoriza- 
tion of Seiberg-Witten curves with random matrix considerations. Another 
interesting open problem would be to generalize the solution to the case of 
non-complete factorization i.e. with at least one monopole staying massive. 
This is of direct relevance to the study of the global structure of A/" = 1 
vacua [32, 33, 34]. 
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